Lifetime of Oil Drops Pressed by Buoyancy Against a Planar Interface. II. Large 

Drops 
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Abstract 

In a previous report [C. Rojas, G. Urbina-Villalba, M. Garcfa-Sucre, Phys. Rev. E 81, 016302 (2010)] it was shown 
that Emulsion Stability Simulations (ESS) are able to reproduce the lifetime of micrometer-size drops of hexadecane pressed 
by buoyancy against a planar water/hexadecane interface. It was confirmed that small drops (r; < 10 pm) stabilized with 
/3— casein behave as non-deformable particles, moving with a combination of Stokes and Taylor tensors as they approach the 
interface. Here, a similar methodology is used to parametrize the potential of interaction of drops of soybean oil stabilized 
with Bovine Serum Albumin (BSA). The potential obtained is then employed to study the lifetime of deformable drops in the 
range 10 /im < < 1000 /^m. It is established that the average lifetime of these drops can be adequately replicated using 
the model of truncated spheres. However, the results depend sensibly on the expressions of the initial distance of deformation 
and the maximum film radius used in the calculations. The set of equations adequate for large drops is not satisfactory for 
medium-size drops (10 /^m < < 100 /^m.) and vice versa. In the case of large particles, the increase of the interfacial area as 
a consequence of the deformation of the drops generates a very large repulsive barrier which opposes coalescence. Nevertheless, 
the buoyancy force prevails. As a consequence, it is the hydrodynamic tensor of the drops which determine the characteristic 
behavior of the lifetime as a function of the particle size. While the average values of the coalescence time of the drops can 
be justified by the mechanism of film thinning, the scattering of the experimental data of large drops cannot be rationalized 
using the methodology previously described. A possible explanation of this phenomenon required elaborate simulations which 
combine deformable drops, capillary waves, repulsive interaction forces, and a time-dependent surfactant adsorption. 

PACS numbers: 82.70.Dd, 82.70.Kj, 82.70.-y, 47.55.D-, 47.57.Bc, 87.14.E-, 07.05.Tp 
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I. INTRODUCTION 

Proteins are polyelectrolytes with a complex structure, 
whose spatial conformations strongly depend on the pH 
of the dispersing medium. They are commonly used as 
stabilizers in food emulsions, although it is often difficult 
to justify the stability of these emulsions in terms of the 
molecular properties of the proteins. The adsorption of 
these molecules to the oil/water interface is usually slow, 
markedly affecting the viscoelastic properties of the in- 
terface and the coalescence time between the drops. 

It is not yet possible to measure the coalescence time 
between two oil drops of micrometer size suspended in 
quiescent media in the absence of external forces. How- 
ever, the diffusion tensor of the drops depends on the ap- 
plied force. In one approximation to the problem, exper- 
imental measurements of the coalescence time between a 
drop and a planar interface are used. Proteins are ex- 
pected to favor the occurrence of tangentially immobile 
interfaces, due to their high molecular weight. Dickinson 
et al. [1] showed that the lifetime of the drops of hex- 
adecane stabilized with /3-casein, K-casein and lysozyme 
decreases with the radius of the drops (l/xm < < 5 ^m) 
when they are pressed by buoyancy against the interface. 
These findings were successfully explained by Basheva 
et al. [2] arguing that when a drop is released in the 
bulk of a liquid, it moves according to Stokes law until it 



approaches the interface to a sufficiently small distance. 
Then it markedly decelerates due to the increase of the 
viscous friction in the remaining gap. If the drop retains 
its spherical shape while approaching the interface, its ve- 
locity can be expressed as a combination of the diffusion 
tensors of Stokes and Taylor, respectively: 
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1 



(1) 



where v^t = F/Gnriri refers to the Stokes law for motion 
of a sphere in an unbounded liquid, and uxa is given by: 



hF 

Gnrirf 



(2) 



In Eq. (2) h is the closest distance of approach between 
two spheres, rj is the dynamic viscosity of the external 
liquid, ri is the radius of the small droplet, and F is the 
external driving force: 



F = -TrrfApg, 



(3) 



where Ap is the density difference, and g is the gravity. 

The time elapsed from the moment in which the drop 
starts to move slowly, until it coalesces with the large 
homophase, is equal to: 
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where v{h) is the velocity of thinning of the hquid film 
between the drop and the interface, h is the closest dis- 
tance of approach between their surfaces, /ijni is the gap 
width at which the thinning begins, and hait, is the min- 
imum distance that can be attained (critical thickness of 
rupture) before the film breaks and coalescence occurs. 

Substituting Eq. (1) into Eq. (4), an inverse depen- 
dence between the lifetime of the drop at the planar in- 
terface, and its radius is obtained: 



T = 
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log 
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(5) 



Recently, our group designed a methodology that 

proved successful in reproducing the experimental data of 
Dickinson et al.. Among other results, it was found that: 
(a) micrometer-size drops behave as non-deformable par- 
ticles which move with a combination of Stokes and Tay- 
lor tensors as they approach the interface; (b) The coa- 
lescence time of the simulations can be fitted with high 
accuracy to the equation: t = 1/{A + B rj), where A = 
-0.054103 s-\ B = 0.075662 (/xms)-i (r^ = 0.9997). 
An approximate power-law dependence indicates that: 



T oc r. 



-1.39 



(r = 0.9949) for this size range; (c) the po- 
tential of interaction has a significant infiuence on shape 
of T VS. Vi. The experimental data can only be accu- 
rately reproduced assuming negligible repulsive barriers; 
and (d) for small drops the uncertainty in the exact form 
of the total potential of interaction between a drop and 
the interface in the presence of proteins mostly resides 
in the steric contribution. This part of the potential de- 
pends on the Flory-Huggins parameter, the interfacial 
area of the protein, and the thickness of the protein layer 
which lies above; the interfacial boundary. 

In order to replicate the lifetime of larger drops 
(10 /um — 1000 /xm) published by Basheva et al. [2], we 
had taken advantage of the dependence of the coalescence 
time on the steric potential. First, we parametrized the 
steric potential of micrometer-size drops by reproducing 
the their experimental lifetime using the tensors of Stokes 
and Taylor. Then, we investigated the dependence of r 
vs. ri for bigger drops employing the parameters of the 
potentials previously selected. 

Unlike small drops, large drops are expected to deform 
as they approach the interface. According to Reynolds 
[3] the velocity of drainage of a thin liquid film between 
two circular sections can be approximated by: 



Re 



(6) 



where rf is the radius of the film. This radius can be 
found from the stress balance of the film at the planar 
surface [2]. It increases with the radius of the drop and 
inversely proportional to the interfacial tension (7): 



Tf — 



7r7 



(7) 



Substituting Eqs. (6) and (7) into Eq. (4) an expression 
for the lifetime of deformable droplets is obtained [2] : 



T = 



47^72 
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(8) 



It is clear from Eq. (8) that in the case of large drops, 
the lifetime increases as a function of a particle radius. 
However, the force usually depends on the particle ra- 
dius, and the critical distance of rupture depends on the 
radius, the interfacial tension, and the force. An estima- 
tion of the coalescence time can be obtained discarding 
the second term of Eq. (8), and using the equation pro- 
posed by Vrij [4, 5] for the critical thickness of rupture. 

In this case, r oc r^^^''' . 

This article is structured as follows: In Sec. II an 
overview of Emulsion Stability Simulations is presented. 
In Sec. Ill the computational details of the simulations 
are given and the parametrization procedure explained. 
In Sec. IV the results are shown and in the final section 
V the conclusion is given. 



II. EMULSION STABILITY SIMULATIONS 

A detailed description of the algorithm of Emulsion 
Stability Simulations can be found in [6 10] and refer- 
ences therein. In ESS the particles move with an equa- 
tion of motion similar to the one of Brownian Dynamic 
Simulations [11]: 

rp,i {t + At) = rp,i (t) + ^^At + R, (9) 

where Tp^j is the position of particle i, Di is the diffu- 
sion constant, is the total force acting on z, kg is the 
Boltzmann constant, T the temperature. At the time 
step, and R a random term which represents the Brown- 
ian motion of the particle. The diffusion constant is equal 
to Di = i'o/corr/corr whcrc Dq = ksT /GiTrjri. The first 
correction term, /iorr, takes into account those factors 

that change the diffusion constant of a particle at infin- 

(2) 

ity dilution. The second correction term /corr takes into 
account the hydrodynamic interactions between the par- 
ticles, caused by the movement of the surrounding liquid 
as the particles diffuse [8]. 

A calculation begins distributing a set of oil drops in 
a cubic box of side length L. It is assumed that the 
molecules of oil mainly determine the van der Waals in- 
teraction between the particles. Instead, the repulsive 
interactions depend on the amount and chemical nature 
of the surfactant molecules adsorbed to the interface of 
the drops. The program has several routines for appor- 
tioning surfactant molecules amongst the drops. In the 
most simple case, the surfactant is distributed evenly and 
instantaneously between the available interfaces. Once 
the surfactant has been allocated, the surface properties 
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of the drops (like charge, interfacial tension, etc.) can 
be computed. Following, the diffusion constant and in- 
teraction forces can be calculated and the drops moved 
according to Eq. (9). At every time step, the program 
checks for the coalescence of drops. In the case of non- 
deformable drops, coalescence occurs whenever the dis- 
tance of separation between the centers of mass of the 



drops, 



is smaller than the sum of their radii. When 



this happens, the former drops disappear and a new big- 
ger drop is created at the center of mass of the colliding 
particles. The radius of the new drop results from the 
conservation of volume. 

The present version of the code can simulate the be- 
havior of systems of non-deformablc or deformable drops. 
In both cases the particles follow the same equation of 
motion (Eq. (9)) but the analytical form of the diffusion 
tensors and the interaction forces change. 

If the mode of deformable droplets is selected, it is as- 
sumed that the deformation occurs independently of the 
energy required for this process. The model of truncated 
spheres is used to simulate the change of shape of the 
drops [12-14]. According to this model, three regions of 
approach can be defined: 

Region I: The distance of separation between the centers 
of mass of the drops, rij, is larger than -|- rj + 
hini, where /ijni stands for the initial distance of 
deformation of the drops. Hence, the drops behave 
as spherical particles. 

Region II: This region covers the range of distances be- 
tween the beginning of the deformation n ^ 0, 
and the attainment of the maximum film radius: 
ff = r-fmax- As soon as the drops enter Region II, 
they change their shape from spheres to truncated 
spheroids. Following the model proposed by Danov 
et al. [13, 14] the closest distance of separation be- 
tween the surfaces of the drops is assumed to be 
constant h = /lini [12] while the radius of the film 
increases: 



^■\/^f ^fmax + \Af ''fmax^ 



hini + \ y I i ' fmax 

^ ^ij ~t~ ^ini^ 



(j^ij ^ini) 



n + ri 



(10) 

(11) 



Region III: The maximum film radius has been attained, 

= ^fmax, and the intervening liquid between the 
drops drains until it reaches a critical distance of 
rupture: 



fmax 



fmax 



(12) 



< /iini +[^Jrf- r-Lax + ^/^l - ^fmax) ' 

= - (V^'-r-La^ + V^f-^w) • (13) 



Accurate estimation of the initial distance of deforma- 
tion /ijni, is very difficult since it results from a balance 
between hydrodynamic and interaction forces. In the 
present calculations the soundness of two diff'erent ex- 
pressions for /lini was compared. The first equation was 
obtained fitting the curves of hini{ri,^) -given in Ref. 
[14]- with a polynomial expression: 



= [1.2932 X 10^ - 8.6475 x 10"^ 

X exp(-r,/1.8222 x 10"*^] 

3.3253 -t- 5.9804 cxp(-7/0.00402) 



3.3253 -I- 5.9804 exp(-10-3 /0.00402) ' 



When the above expression is used to estimate hi, 
maximum film radius is approximated by: 



rfn 



(14) 



i, the 



(15) 



As an alternative formula for hmi, we considered a frac- 
tion [fh) of the expression obtained using the lubrication 
approximation in the presence of the buoyancy force [15]: 



hi. 



f, 



2rfApg 
37 



(16) 



where fh is a real number between zero and one. When- 
ever Eq. (16) is used, the maximum film radius is esti- 
mated as a fraction {fr) of the value predicted by Eq. 
(7): 



r-fn 



l2gAp 
37 ■ 



(17) 



When either an arbitrary estimation of /lini is used as 
an input of the simulation, or an alternative equation like 
Eq. (16) is used for its evaluation, Eq. (11) cannot be 
employed to estimate Vf. For this purpose an alternative 
expression is implemented: 



{nj - hini f - (r| - rf ) 

2 (j'ij ^ini) 



n 2 



(18) 



Eq. (18) can be deduced assuming that the film formed 
between two drops has a uniform thickness. As a conse- 
quence, the radius of the deformation in each drop is the 
same, and: 



(19) 
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where dfe is the distance between the center of the trun- 
cated sphere k and the center of its planar interface. No- 
tice that when Eqs. (18) and (19) are used, the surface- 
to-surface separation between the drops (h) changes dur- 
ing the growth of the film radius (Region II). 

In regard to /icrit, the expression published by Sche- 
ludko and others was used [16-19]: 



where Adit — rf/lO and Ah is the Hamaker constant. 

Even in the mode of dcformable drops, the particles 
behave as spheres if > r.^ -I- rj + hmi- This means 
that the potential of interaction and the diffusion con- 
stant correspond to the ones of spherical particles within 
Region I. At h = /lini, the code calculates the dimensions 
of truncated spheres which are compatible with the ac- 
tual distance of separation between the centers of mass 
of the spherical drops {rij < ri + rj + /lini)- In this case, 
the expressions of the potentials corresponding to two 
truncated spheres are employed to move the drops (see 
Table I and Ref. [14]). Those potentials are expressed in 
terms of the thickness of the film and its radius. The use 
of three regions of approach allows to develop equations 
which relate r^j to those variables. As a result, the po- 
tentials can be differentiated algebraically, and the force 
between the particles can be calculated (see Refs. [9, 10] 
for details). 

The total potential of interaction between a non- 
deformable drop and the interface, is assumed to be 
composed of three contributions: van der Waals, electro- 
static, and steric. In the case of deformable drops, two 
new potentials of interaction appear during the evolution 
of the film (Region II). They take into account: (a) the 
surface deformation energy (extensional or dilational en- 
ergy) coming from the increase of interfacial area as the 
spherical drops turn into truncated spheres, and (b) the 
bending elasticity of the surfactant monolayer adhered 
to the interface of the drops [12]. These two potentials 
change with the interparticle distance during the forma- 
tion of the film (Region II), but reach a constant value 
after a maximum film radius has been attained. Hence, 
they do not contribute to the value of the force within 
Region III. 

In order to simulate the movement of the drops, three 
diffusion tensors were used (Table II). In the case of non- 
deformable drops the expressions of Stokes and Taylor 
were implemented. In the case of deformable drops, the 
expression of Danov et al. [13] for deformable drops was 
used. In the last case, the value of the parameter 
was fixed to 1.0, in order to simulate the behavior of 
tangentially immobile interfaces. 



III. COMPUTATIONAL DETAILS 

The parameters employed in the simulations are shown 
in Table III. They correspond to a soybean in water 
emulsion stabilized with BSA. The protein concentration, 
pH and ionic strength of the aqueous solution correspond 
to the the experimental measurements of Basheva et al. 
(4 X 10"" wt% protein, pll=6.4, ionic strength 0.15 M 
[2] ) . The Hamaker constant was approximated using the 
expression of Lifshitz for the case of two identical slabs 
of oil separated by water [26] . In the absence of more ac- 
curate data, the refractive index, dielectric permittivity, 
and the main absorption frequency in the UV of olive oil 
were used for this purpose [27, 28]. The electric charge 
of the BSA molecule was calculated by reproducing the 
value of the C, potential (C=— 47 mV) of a 0.5 /im soybean 
particle covered by BSA (Kong et al. [29]). The inter- 
facial area of the protein r^J,^^ was also calculated from 
the data reported on Ref. [29] . The width of the protein 
layer was obtained from Ref. [30] (see also [31, 32]). 

Following the methodology employed in Ref. [10], the 
fiat oil/water (0/W) interface was represented by a very 
large drop of oil fixed in the space (Fig. 1). A preliminary 
set of calculations was necessary in order to establish the 
radius of the large drop and the time step of the sim- 
ulations. For that purpose, the Flory-Huggins solvency 
parameter x was temporarily set to 0.4 (for globular pro- 
teins it usually varies in the range 0.3-0.5 [33]). Radii of 
500 /im, 5000 /im, and 10, 000 /im were tested for the large 
drop, and scaled time steps (At* = AtDo/rf) of 0.005, 
0.01, 0.05, 0.1 and 1.0 were studied. The volume fraction 
of protein around the drops (</>) was used as an input of 
the simulation (0.0056 < 0i < 0.0670) for several initial 
distances of separation between the interface (large drop) 
and the approaching droplets (30 nm < d < 500 /im). 
The last two variables were systematically changed until 
the order of magnitude of the experimental coalescence 
time was reproduced for drops smaller than < 15 /im. 
It was confirmed [10] that the repulsive potential barriers 
between the drops had to decrease with the radii of the 
drops in order to obtain a set of coalescence times that 
diminish as a function of the particle radius. 




~ ^ Oil droplet 
of radius r 

Real System Model 

FIG. 1: Model of the drop/interface system employed in the 
simulations. Here, r is the radius of the small droplet, and R 
the radius of the large drop resembling the interface. 
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TABLE I: Potentials used in the simulations. In these equations, r, is the radius of the small droplet, Vj is the radius of the 

large drop and h is the minimum distance between their surfaces. For the van der Waals (vdW) potential: Ah is the Hamaker 

constant, x = y = I = h + rt + y^r'f — r'f, L = rj + — rf, d = \Jh? -\- Ar'f, h and rt arc the thickness and radius 

of the film, respectively. For the electrostatic potential (elect): = ^Jfe^T ^'^'^ charge number, e is the dielectric 

permittivity of the medium, Cd is the electrolyte concentration, e is the electron charge, fcsT is the thermal energy, ^'si and 
are the surface potentials for the small and large drops, respectively. For the steric potential (st): is the molar volume of 
the solvent, x is the Flory-Huggins solvency parameter of the protein, 4>j and 4>i are the average volume fraction of the protein 

around each sphere, 4>i = — ^ r ^' '"^I'g — jt, with F the number of molecules per unit area, pp the density of the protein, Mp 

the molecular weight of the protein and 5 the width of the protein layer. Volumes Va, Vb, and Vc depend on h. Their explicit 
geometrical expressions can be seen in Ref [20]. For the dilatational (extensional) potential (dil): 70 is the interfacial tension 
and ra = ^rip^- For the bending potential (bend): |Bo| = 1.6 x 10~^^ N [21, 22] is the interfacial bending moment. 



Kdw = -4f- [ + '2^ \ ^ +2ln( f+^'l+l )] . [23] 

12 ^x'^+xy+x x-^+xy+x+y \x-^+xy+x+y J j L J 



Spherical 



V,i,,t = 2^CeifesTtanh 



tarih 



3Vi 



(i-x)(5-tr[ 



h+ri+rj)\ ' 

S <h<2S. 
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[14] 

[20, 24] 



v.. = ^ (I - X) [i^^r a - Va) + (0.)' (S - + 20.^i (^)] , [20] 

0<h<6. 



VvdW - - — l(L+h) + h{l+L) + ^"1 [l(h+L)\ + W 

2rf _ i-r,-(L-rj) 2r| _ 2(L-r^)-fe 2rj L^(i-fc) 

L hi 2l-2ri-h hi 2l-2ri 2h ' hl(l + L)(L + h) 

2r| l^+rf 2r|d 

h(2l-2ri-h) (l+L)(l+L-2rj) + (2i-2ri -h) [(h+L){h+I,-2r-j- )-(i-h)(i-2ri -Jt)J . 
4r?{I-h) 

(i+i)(i+-£'-2rj)[(/i+L)(h+L-2rj)-((-ft,)((-2ri-;i)] 




Deformable 



Select = 2^C.ifeBTtanh (^) tanh (3^) x e''^'^ [rj + 

(i-x)(5-tr[^+2<5+t-toii], 



3Vi > 



[14] 
[20] 



5 <h<25. 



= ^ (I - X) [{hf - + (I - -a) + 20.^, (i^)] , 

< /i< 5. 



[20] 
[14] 



Vbend = -2^, (rf/r„)2 « 1. [12] 



Once the values oi d {d = 35 /xm). At* (see Table IV), 
and rj (rj = 5000 /um) were fixed, the program was al- 'ir^TM 

lowed to re-calculate ^, using the following relation: = ^, (21) 

PpNa [{n + Sy - rf] 

where pp is the density of the protein, Mp its molecular 
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TABLE 11: Tensors used in the simulations. Here Vi is the radius of the small droplet, rj is the radius of the large drop, 
— r^+r- ' ^ ^^'^ ^^'^ thickness and radius of the film, Do = q^^. , V is the dynamic viscosity of the continuous phase 
and es have values between 0.001 and 1. 



Geometry 


Tensor 


Stokes Immobile: A sphere in an 


Dst = 


-Do. 






unbounded liquid [2]. 










Taylor Immobile: Two spheres of radii 


I>Ta = 


= 4Do 


ri 




T'i.rj, and h <C 'ri, rj [25]. 










Two deformed drops [13]. 


Dud = 


4h 

n 


(- 


2 4 \ -1 

nh rfh?) °' 



TABLE III: Parameters of the simulations. The molecular 
properties correspond to one BSA protein. 

Hamaker constant [26] 2.72 x lO'^^ J 

Ionic strength [2] 0.15 mol/1 

Surf, concentration [2] 6.04x10"* mol/1 

Molecular mass [34] 66.267 kDa 

Molar volume [34] 0.048574 m^/mol 

Electric charge — 6.135e 

Area 11007.81 

5 [30] 1.1 nm 

X 0.49993 

7 [2] 15 mN/m 



weight, and 5 the width of the protein layer around the 
drops. The new values of 4>i predicted by the Eq. (21) 
were considerably larger than the ones formerly tested 
(for example, = 0.662 for = 6.9 y^m). Hence, the 
value of X was now changed until the order of magnitude 
of the coalescence time coincided again with the one of 
the experiments. Prom this procedure, a value of x = 
0.49993 was obtained. 

Finally, the values of At,rj,(j)i, and x optimized for 
non-deformable particles were then used without further 
modification for the calculations of dcformable droplets. 
However, when Eqs. (14) and (15) were employed to 
evaluate /lini and rfmax, deformablc drops required values 
of d of the order of nanometers in order to reproduce the 
order of magnitude of the experimental coalescence time. 
Instead, non-deformable drops needed initial separations 
of a few microns. 

Figures 2 and 3 show the potentials of interaction ob- 
tained for non-deformable and dcformable drops, respec- 
tively. The potential of dcformable drops shown in Fig. 
3 was calculated by using Eqs. (16) and (17) for hini and 
^fmax- Notice that the potentials shown in these figures 
are plotted as a function of h to highlight the location of 



TABLE IV: Time step of the simulations. 





r (pm) 


At* 


At (s) 




1-7.5 


1 


2.0 X 10 


-'^ - 8.6 X 10" 


-4 


Spherical 


10-20 


10-' 


2.0 X 10 


- 1.6 X 10" 


-4 




25-40 


10-2 


3.2 X 10 


- 1.3 X 10" 


-5 




80 - 200 


10-^^ 


1.0 X 10 


- 1.6 X 10" 


-3 




25-90 


10-2 


3.2 X 10 


- 1.5 X 10" 


-4 


Deformable 


100 - 375 


10-^^ 


2.0 X 10 


- 1.1 X 10" 


-4 




400 - 1000 


10"^ 


1.3 X 10 


- 2.0 X 10" 


-7 



the repulsive barrier (whenever it occurs). As a result, 
the width of Region 11 cannot be appreciated. In the 
potentials shown in Fig. 3, the tiny range of distances 
corresponding to the abrupt jump in potential, marks the 
growth of the film radius (Region II). 

In the case of non-deformable particles, the total po- 
tential of interaction is always attractive, increasing in 
absolute magnitude as a function of the particle radius. 
Instead, the potential of deformable droplets shows a very 
peculiar behavior. For a drop with radius < 100 /im, 
it decreases with the increase of the radius. However, for 
ri > 100 /Ltm a repulsive barrier develops. This barrier 
increases with the particle radius. 

A close look at the partial contributions of the poten- 
tial of deformable drops indicates that the repulsive bar- 
riers illustrated in Fig. 3 are caused by the cxtcnsional 
contribution. At small radii, all the potential contribu- 
tions occur at a similar distance of separation. Within 
a few nanometers, the repulsive contributions decrease 
in the following order: electrostatic > steric > bend- 
ing > extensional. However the van dor Waals potential 
prevails for n < 100 /xm and the total potential of in- 
teraction is attractive at all distances. As the size of the 
drops increases, /ijni also grows. Thus, Region II progres- 
sively moves toward longer interparticle distances. The 
range of action of these potentials enlarges (Region III) 
but the force is exerted within Region II. Thus, the range 
of action of the extensional and the bending force sepa- 
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rates considerably from the one of the electrostatic and 
the steric contributions. Moreover, the cxtensional con- 
tribution gradually surpasses the bending contribution. 
At Tj = 100 /xm, the electrostatic potential shows a peak 
around 9 nm, while the extensional contribution reaches 
its maximum around 26 nm. For — 200 /im the ex- 
tensional barrier already occurs a.t h = 208 nm. Since 
the absolute magnitude of the van der Waals potential 
increases with both and rf, the electrostatic and the 
steric repulsions are always suppressed, and the total po- 
tential is attractive at short separations. Yet, the ex- 
tensional and bending repulsions grow within Region II, 
exceeding the van der Waals attraction at longer separa- 
tions (see Fig. 3). From that point on, a repulsive barrier 
can be observed in the total potential of interaction. 




.•a -600 - 




h (nm) 

FIG. 2: Total potential of interaction between a large spher- 
ical drop of 5000 /xm and a small spherical droplet of micron- 
size. Dotted line: r = 2 /um; dot-dashed line: r = 4 /xm; 
dashed line: r = 7.5 /xm; dot double-dashed line: r = 10 /um; 
dash double-dotted line: r = 25 //m; solid line: r = 100 /xm. 

The present simiilations consist in the evaluation of 
the average coalescence time between a set of drops of 
different radii and the interface (sec Fig. 1). The aver- 
age coalescence times were computed from 300 random 
walks. For preliminary simulations a selected number of 
particle radii were explored. For the final calculations 
an extensive number of radii were computed. From 1 /um 
to 4.75 /Lim the radius was changed using increments of 
0.25 /im. From 10 /im to 50 /im increments of 5 /im were 
utilized. From 50 /xm to 100 jim, 10 /xm-increments were 
employed. Finally, for the range 125 /xm to 1000 ytxm, in- 
crements of 25 /im were used. For each particle radius the 
average coalescence time was first evaluated for spherical 
particles, and then repeated for deformable drops. 

The experimental data in the range 1.48 /xm to 234 /xm 
was kindly provided by Dr. T. D. Gurkov. The data 
corresponding to larger particle radii, was obtained di- 
rectly from Fig. 2 of Ref. [2] using the program Engauge 
Digitizer. 




h (nm) 

FIG. 3: Total potential of interaction between a large de- 
formable drop of 5000 /im and a small deformable drop of in- 
termediate size (/lini and rfmax were calculated with Eqs. (16) 
and (17), respectively). Dash double-dotted line: r — 25 /xm; 
dashed hne: r = 50 /xm; solid line: r = 100 /xm; dot-dashed 
line: r = llO/im; dotted line: r = 120 /ttm; dot double-dashed 
line: r = 130 /xm. 



IV. RESULTS AND DISCUSSION 

A. Average Value of the Lifetime of the Drops as 
a Function of the Particle Radius 

The apparatus used in the experiments published by 
Basheva et al. [2] essentially consisted on a glass vessel 
with a planar oil/water, and a thin glass capillary sit- 
uated underneath, 0.5 cm below the surface. Different 
procedures were applied for the making of small drops 
(micrometer size up to rj < 100/xm), medium size drops 
(100 /im < Ti < 500 /xm) and large drops (r^ > 500 /xm). 
Micrometer size drops were produced by preparing an 
emulsion with a rotating blade homogenizer. Thus, they 
were introduced in the experimental cell by means of a 
glass capillary and a syringe. Medium size droplets were 
formed squeezing out a drop of oil from the capillary, and 
then, sucking back the oil in a sudden and fast manner. 
This produced a turbulent flow of oil and water that en- 
ter the capillary, favoring the formation of an emulsion 
in situ. That emulsion was pushed out subsequently. For 
large drops a single drop of oil was created at the tip of 
the capillary and then it was blown out by a certain vol- 
ume of solution with the help of an additional glass tube. 
A CCD camera mounted on a microscope was used for 
recording the process. A timer allowed determining the 
droplet lifetime, measured from the moment in which the 
oil droplet begins to move slowly until it coalesces with 
the bulk oil phase. 

Unlike the above definition given by Basheva et al. [2] , 
we will use the terms "medium-size" or "intermediate- 
size" to refer to those drops in the range (10 /xm < < 
500 /xm). 

Fig. 4 illustrates the results of the calculations for 
non-deformable micrometer drops which showed the best 
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agreement with the experimental data. The hnes in the 
figure are only a guide to the eye. As previously dis- 
cussed, the data of micrometric drops was used to adjust 
the parameters of the steric potential. Notice that the 
data of Bahseva et al. [2] has a considerable dispersion 
at all ranges of particle sizes (see below). Due to this fact, 
it is not possible to reproduce the points corresponding 
to drops of micrometer size with the same accuracy pre- 
viously achieved for the experimental data of Dickinson 
et al. [1]. The theoretical curve was produced using a 
fixed distance between the interface and the initial posi- 
tion of the moving drops {d = 35 /im, solid line in Fig. 4). 
It is observed that the lifetime of small drops decreases 
as a function of the oarticle size. 
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FIG. 4: Average coalescence time vs. droplet radius for spher- 
ical non-deformable drops. Crosses: experimental data [2]; 
solid line: constant initial distance of separation between the 
drops and the interface (d = 35 //m); dashed line: variable 
initial distance of separation between the drops and the in- 
terface, d('r). See text for details. 

The use of a distinct initial distance of approach d{ri) 
for drops of different sizes (dashed line), favors a better 
quantitative agreement between some of the points of the 
experimental data and the simulations (Fig. 4). How- 
ever, it increases the difficulty of the calculations con- 
siderably, introducing an undesirable additional variable 
that needs to be adjusted. This point is illustrated in Fig. 
4 using a limited number of simulations. In order to im- 
itate the effect of a variable d(ri), Eq. (5) was equalized 
to the empirical expression of the lifetime of the drops 
found by Basheva et al. [2]: t = O.I45/r, [2]. The nu- 
merical solution of the resulting algebraic equation (using 
Eq.(23) for ft-crit) yielded a value /lini for each particle ra- 
dius. The values of /lini changed from 17 to 1200 /xm for 
(1 jim < ri < 100 iJ,m). Following, we approximated d{ri) 
by /lini and ran simulations of non-deformable drops. Ob- 
serve that d represents the distance below the interface 
at which the small drops are released, while hini stands 
for the separation at which the analytic form of the ten- 
sor changes from Stokes to Taylor. The error bars of 
the theoretical points were approximated by the stan- 



dard deviation of 300 random walks. In order to optimize 
the clarity of the figures, the error bars are only shown 
in Fig. 4. The uncertainty of the simulations increases 
with the decrease of the size of the drops, evidencing the 
effect of the Brownian motion on the movement of the 
drops. Notice also that the error bars are already small 
for medium size drops. This indicates that the scattering 
of the experimental points corresponding to large drops 
is not caused by their Brownian motion. 

As can be appreciated from the inset of Fig. 4, the 
curve corresponding to spherical drops also decreases 
with the particle radius for medium-size drops. Such 
behavior is essentially the consequence of two factors. 
First, the attraction between the drops and the interface 
increases with the size of the drops. Second, the Tay- 
lor tensor (see Eq. (2)) is inversely proportional to the 
square of the particle radius but directly proportional to 
the buoyancy force. Hence it increases with the radii of 
the drops, favoring a hyperbolic decrease of the lifetime 
of small drops as a function of their size. Consequently, 
the increase of the coalescence time observed for large 
drops as a function of their radii, cannot be explained 
assuming non-deformable drops. 
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FIG. 5: Average coalescence time vs droplet radius for de- 

formablc droplets. Crosses: experimental data [2]; solid line: 
d = 35 pim; dashed line: non-deformable drops [d = 35 /im); 
dot double-dashed line: d, = 30 nm; dot-dashed line: d = 
50 nm; dotted line: d = 60 nm; dash double-dotted line: 
d = 100 nm. 

Fig. 5 illustrates similar calculations for deformable 
drops. In this case, the effect of the initial distance 
between the drop and the interface is illustrated using 
distinct values of d. In all calculations of deformable 
droplets shown in this paper the tensor of Danov et al. 
[14] was used. However, different equations were used 
to approximate the initial distance of deformation and 
the value of the maximum film radius. In the simula- 
tions of Fig. 5, Eqs. (14) and (15) were employed to 
estimate /lini and rfmax- It is observed that initial dis- 
tances of approach (rf) of the order of tens of nanometers 
(20 nm— 100 nm) are sufficient to attain the order of mag- 
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nitude of the experimental data. The change in the value 
of d with respect to the one of smaller non-deformable 
drops is due to the dynamics of thinning of the aqueous 
film between the drop and the interface. The thinning 
of a plane parallel film is very slow. Hence, a spherical 
drop needs to diffuse through a longer distance in order 
to spend a similar period of time before coalescing. 
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FIG. 6: Average coalescence time vs droplet radius for de- 
formable drops (/iini was calculated using Eq. (14)). Crosses: 
experimental data [2]; solid line: deformable drops (d = 
35 fim.) 

Unfortunately, the region between 10 /xm < < 30 /um 
does not contain sufficient experimental points to vali- 
date which type of drop is more likely to occur. For drops 
of intermediate sizes the number of experimental points 
that can be justified using simulations of deformable 
drops is considerably larger. However, the curves of de- 
formable drops do not overlap smoothly with ones of 
spherical droplets at small particle radii. Moreover, they 
show an abnormal decrease close to the hyperbolic rise 
of the coalescence time, exhibited by the experimental 
points and the simulations of non-deformable particles. 

Changes in d from 30 nm and 100 nm improve the 
agreement between the experimental data and the cal- 
culations. Such a variation appears feasible due to the 
procedure employed for making the drops of micrometric 
size. An emulsion is first prepared, and their drops are 
released from a capillary below the interface. In fact, ac- 
cording to Basheva et al. [2] the experimental procedure 
allows to measure the size of the drops with an accuracy 
of 0.5 fiui. Therefore, it is unlikely that changes in the ve- 
locity of approach corresponding to separation distances 
(d) lower than 500 nm can be determined. 

When the predictions of the simulations of spherical 
and deformable droplets are extended to larger particle 
sizes (10 /im to 500 /im), the curves of r vs. mono- 
tonically decrease (Figs. 4, 5 and 6). In the case of 
deformable drops (inset of Fig. 5 and Fig. 6), this is 
caused by the expressions of h-mi and rf,„ax used in the 
simulations. When the radii of the particle increases be- 
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FIG. 7: Total potential of interaction between a large de- 
formable drop of 5000 /jm and large deformable drop of a few 
hundred microns (/lini was calculated using Eq. (14)). Dot- 
ted line: r = 25 /xm; dash double-dotted line: r = 100 /xm; 
dot-dashed line: r = 300 /xm; dashed line: r = 500 /tm; solid 
line: r = 800 /tm. 



yond 100 /xm, the variation of /iini(?"i)7) vs. predicted 
by Eq. (14) approach an asymptotic limit of a few tens 
of nanometers. As a consequence, the increase of rfmax 
as a function of is also dampened, and the potential 
of interaction decreases with the increase of the parti- 
cle radius (Fig. 7). However, it is not the increase of 
the van der Waals attraction which causes the decrease 
of the coalescence time as a function of the particle ra- 
dius. Figure 8 shows the results of a limited number of 
simulations of deformable drops in which no interaction 
forces were used. For this case, Eqs. (16) and (17) were 
employed. Neither the extcnsional nor the bending po- 
tentials were included in the simulations. In this case, the 
lifetime of the drops increases with the size of the drops. 
Moreover, if the same simulations arc run including the 
van der Waals potential (solid black circles in Fig. 8) , an 
analogous result is obtained. This illustrates the effect of 
the hydrodynamic tensor on the movement of the drops. 
It also evidences the influence of the expressions of hmi 
and Tfmax in the calculations (see below). 

Fig. 9 shows the results of similar simulations in which 
the total potential of interaction was included. This po- 
tential increases with the particle radius as a consequence 
of the extensional contribution (Fig. 3). As in the pre- 
vious calculations, Eqs. (16) and (17) were used as ap- 
proximations for /lini and rtmax- The coalescence time 
increases pronouncedly with the radii of the drops. Sur- 
prisingly, the increase of the lifetime of the drops as a 
function of the particle radius is not caused by the de- 
velopment of the repulsive potential barrier. This can be 
realized comparing the dashed curve and the black cir- 
cles illustrated in Fig. 8: the lifetimes predicted by the 
simulations in the presence or absence of repulsive forces 
are similar. 

Table V shows a comparison between the largest mag- 
nitude of the repulsive force (caused by the extensional 
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FIG. 8: Dependence of the lifetime on the potential of in- 
teraction between the drops and the interface. Crosses: ex- 
perimental data [2]; solid line: the movement of the drops 
is determined by their thermal interaction with the solvent, 
and the buoyancy force. In these sirrmlatioiis the potential 
of interaction was coinplctcly suppressed, and the tensor of 
Danov et al. [13] was used, /lini was calculated by means of 
Eq. (16) with fr = 7/20; solid black circles: similar calcular- 
tions but including the van der Waals potential; dashed line: 
similar simulations in which the total potential of interaction 
was included. 
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FIG. 9: Average coalescence time vs droplet radius for de- 
formablc drops. Crosses: experimental data [2]; solid line: 
d = 35 ^nn (Eq. (16) with /h = 1 and Eq. (17) with fr = I 
were used as approximations for /lini and rfmax); dashed line: 
d = 10 /xm with Vf = 0.05 n, and /lini = 10~* m 



deformation of the drops) and the value of the buoyancy 
force. It is observed that in the present case the buoy- 
ancy force subdues the effect of the interaction forces. 
Consequently, our results coincide with the previous pre- 
dictions of Ivanov and Kralchevsky [35] for the case of 
zero disjoining pressure: the lifetime of large drops in- 
creases as a function of their radii (see Fig. 3 in Ref. 
[35]). 

It is important to realize that the extensional potential 
only acts within Region II. The width of this region de- 
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pends on the size of the drops, changing between 30 nm 
and 1.92 /xm in the range 250 ^m < < 1000 /xm. There- 
fore a limited number of additional simulations were run 
decreasing the time step in one and two orders of mag- 
nitude, to improve the sampling of the repulsive barrier. 
However, similar lifetimes were obtained. This fact val- 
idated the magnitude of the time steps used in the sim- 
ulations. Furthermore, it verifies that the extensional 
contribution of the potential was unable to overcome the 
effect of the buoyancy force. 

We also tested several alternative expressions for h\^^ 
and rfmax- The dashed line of Fig. 9 illustrates the re- 
sult of using /lini as an input of the simulations (/lini = 
lO^^m) along with a slight modification of Eq. (15): 
rfmax = frVi'ihini {.fr = 0.055). Somc of these expres- 
sions favor an increase of r as a function of but produce 
a convex curve contrary to the experimental results. It 
is possible that equations of the form: hini oc (rj)", and 
rfmax (r^)™ might bo valuable in the future, but we did 
not make an extensive study of these expressions. In any 
event, it is evident that the predictions of the simulations 
for large drops depend sensibly on the analytic form of 

hini and rfmax- 

In order to reproduce the experimental value of t for 
ri = 1000 /im, the fraction fr was systematically changed 
while keeping fi^ — 1 in Eq. (16). Table VI shows the 
coalescence times obtained for different values of fr- The 
results of this table suggest that a value of fr = 7/20 
might be adequate to reproduce the experimental data 
of large drops. Fig. 10 shows the results of the sim- 
ulations for the complete range of rj using fr = 7/20 
and a constant value of = 35 /zm. This is the value 
of d previously employed in the calculation of spheri- 
cal drops of micrometric size. For particle radii in the 
range 1 — 10 iim we already showed that the simulations 
of non-deformable drops reasonably reproduce the exper- 
imental data. It is now observed that in the range 100 i^m 
< Ti < 1000 fjbm the experimental data is also recreated 
using deformable drops if Eq. (16) is used to approxi- 
mate the initial distance of deformation and Eq. (17) is 
employed to calculate rfmax with the factor /,. = 7/20. 

Despite the above results the extension of the simula- 
tions of deformable droplets to the small range of particle 
sizes, fails to reproduce the experimental data (Fig. 10). 
Since the tensor of Danov et al. tends to the expression 
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FIG. 10: Average coalescence time vs droplet radius for spher- 
ical and deformable drops initially separated by a distance 
of d = 35 frni. Crosses: experimental data [2]; dashed line: 
simulation with spherical drops; solid line: simulation with 
deformable drops. 



potential. At the other extreme of particle sizes, the 
experimental data of micron-size ckops < 10 fiia can 
only be reproduced using attractive potentials and non- 
deformable drops. For these drops the magnitude of the 
attractive potential increases with the size, reinforcing 
the effect of the hydrodynamic tensors (Stokes and Tay- 
lor). Moreover, the effect of the Brownian motion is sig- 
nificant as evidenced by the error bars of the simulations. 
For medium-size particles, the errors bars are very small 
and the outcome of the simulations depends on the inter- 
play between the buoyancy force and the potential of in- 
teraction. However, as the particle radius decreases, the 
deformation of the drop also diminishes. Consequently, 
the repulsive cxtensional potential decreases significantly 
along with the buoyancy force. Hence, the effect of the 
buoyancy force still predominates while the deformation 
of the drops decreases and the hydrodynamic tensor ap- 
proaches the form of the Taylor tensor. 

An empirical fit of the computed lifetime of droplets 
for 300 <ri< 1000 iJ.m, leads to: 



of Taylor for small particle radius, r increases when r, 

decreases. However, this augment occurs at considerably 
larger particle radii than those experimentally observed. 
Moreover, the expression of ft-ini calculated by means of 
Eq. (16), becomes smaller than hait for rj < 25 fim. 
This unphysical result is partially due to the use of Eq. 
(19) to approximate hcru- This equation was deduced 
assuming that only van der Waals forces occur in the 
system. In any event, it is not possible to do simulations 
for rj < 25 jum using Eq. (16) to approximate /lini- 
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FIG. 11: Average coalescence time vs. droplet radius for 

spherical and deformable drops. Crosses: experimental data 
[2]; solid line: simulations of spherical and deformable drops. 
The picture illustrates the results of three different method- 
ologies of calculation for three distinct regions of particle sizes. 
The arrows indicate the overlap between these regions. 

In general, as the particle radius decreases, the effect of 
the buoyancy force is less relevant. The infiuence of the 

buoyancy force on the lifetime of large drops is decisive 
despite the considerable magnitude of the extensional 



T = 3.3494 X 10-^rf (r^ = 0.9985). (22) 

It is remarkable that the value of the exponent in Eq. 
(22) is fairly close to 25/7 = 3.57 suggested by Basheva 
et al. [2]. The theoretical estimation results from substi- 
tuting the critical thickness of rupture suggested by Vrij 
[4, 5, 36]: 



0.268 



^H' f 



6.5Fa 



1/7 



(23) 



into Eq. (8): 



T = 4.088 r?a-«/^A-^/' {^pgf rfl\ (24) 



TABLE VI: Choice of the parameter 
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110.37 


57.10 
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24.42 



As it was explained at the beginning of this section, 
the experimental procedure employed for the production 
of drops of different sizes was distinct for drops of micro- 
metric, medium size, and large size. Fig. 11 illustrates 
the lifetime of the simulations which render the closest 
agreement with the experimental data. For < 2 /im 
spherical drops were calculated using d = 35 /xm. In the 
range 2.5 /zm < < 200 /xm deformable drops were sim- 
ulated using Eq. (14) for ft-ini and d = 50 nm. Finally, 
for the range 225 /xm < Vi < 1000 jim the experimental 
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data was reproduced considering dcformable drops with 
an initial distance of approach of d = 35 fiia {fr = 7/20). 
It is evident that Eq. (14) is more rehable than Eq. (16) 
in order to approximate the value of /lini in the range 
15/Ltm < n < 100 /xm. Since the difference between the 
experimental data and the theoretical prediction is ba- 
sically caused by the method of evaluation of /lini and 
''fmax, it is also reasonable to conclude that the drops be- 
have as dcformable particles over the intermediate range 
of particle sizes. 

For completeness, the change of h vs. t between the 
drops and the interface was calculated for different par- 
ticle radii (Fig. 12). For drops larger than 25 /im, the 
thinning of the film is monotonous and the effect of the 
Brownian movement small. Using the data of h vs. t, the 
velocity of thinning, —dh/dt was evaluated numerically 
at t = 1 s for a selected set of particle radius belonging 
to the region of large drops (r^ > 500 /xm). 



usually involve microscopic distances between film radius 
of macroscopic size (sec Fig. 2 in Ref. [39]). Such ex- 
periments indicate that there is a simple proportionality 
between —dh/dt and h for a wide range of film thickness 
(25 nm to 100 nm) [18]. According to our results: 



dt 



0.9997. 



(26) 



The exponent of Eq. (26) appears to be reasonable 
considering the limitations of the simulations (constant 
d, plane parallel film of uniform thickness, etc.). This 
proportionality between the velocity of thinning and its 
width can be partially justified analytically. If the Brow- 
nian motion of the large drops is assumed to be negligible, 
their velocity is only the product of the diffusion tensor 
and the buoyancy force: 



V = D 



Danov 




FIG. 12: Variation of film thickness with time. Solid line: r = 
25 fim; dashed lino: r = 50 jim; dot-dashed line: r = 100 jim; 
dotted line: r = 300 /im. 

Several models arc available in the literature for the 
calculation of rate of drainage —dh/dt between a drop 
and the interface in the absence of surfactant molecules 
[37]. They reproduce the cases in which: (a) the sphere is 
non-deformable but the interface deforms; (b) the sphere 
and the interface form a plane parallel film; (c) the drop 
and the interface deform forming a curved film; and (d) 
the drop forms a dimple and the interface does not de- 
form [37]. The first three models (a)— (c) predict that: 



dh 
' dt 



ch^ 



(25) 



with a different coefficient c. Despite their limitations, 

the referred formalisms assume the occurrence of tangen- 
tially immobile interfaces as it happens in the presence 
of surfactants molecules. The dependence of —dh/dt on 
h comes from the use of the Stephan-Reynods equation 
[38]. Most experiments that measure /icriti and —dh/dt 



F 



(27) 



Using the tensor of Danov et al. [14] and Eq. (17) to 
relate the film radius with the particle radius gives: 



dh _ cih^ 

dt h'^ + C2h + cs' 



(28) 



The coefficients of Eq. (28) depend on: Ap, g, 7, /,., 
and rj. The relative magnitude of the coefficients deter- 
mine the power dependence of the thinning velocity on h. 
Fovn = 100 /im: ci = 67.18 s^^, C2 = 1.69x lO'i'^ m, and 
C3 = 4.11 X 10~*m^. Likewise, for rj = 1000 /xm: Ci = 
671.82 s-\ C2 = 1.69 X 10-11 m, and C3 = 4.11 x IQ-^m^. 
It is clear then that for the parameters corresponding to 
the experiments of Basheva et al. [2] , the last two terms 
of the denominator of Eq. (28) are negligible and as a 
result: —dh/dt oc h. 



B. Dispersion of Experimental Lifetimes. Possi- 
ble Role of Interaction Forces and Adsorption Times 

It is clear from Figs. 6 and 11 that the experimental 
lifetimes show a considerable scattering at large parti- 
cle radii. Gurkov et al. [25] had pointed out that this 
phenomenon is caused by the statistical nature of the 
coalescence process. The rupture of thin liquid films 
might occur through the spontaneous growth of fluctu- 
ation waves on the two opposing 0/W interfaces of the 
film, or through the formation of holes [40, 41] . In either 
case the rupture of the film is a stochastic process and 
must be treated accordingly. Gurkov et al. [25] studied 
the variation of the number of drops, N(t), which sin- 
vive coalescence until time t (1000 /tm < fj < 3, 500 /tm). 
The relative change {w = 1/N x dN/dt) is related to the 
probability for "drop burst per unit time" . Since w is 
time dependent, it can be expanded in power series: 
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= Wo{t - to) + Wl{t - tof + 



(29) 



where to is the time when coalescence is first observed 
in the system. Surprisingly, use of Eq. (29) produces 



a monotonous parabolic curve of — log vs. t — to, 
where iV^ot is the total number of drops studied (see Fig. 
6inRef. [25]). 

According Ghosh and Juvekar [42], the scattering of 
the rest time of the drops is caused by differences in the 
surface excess of the surfactant corresponding to each 
drop. Those differences promote distinct magnitudes 
of the repulsive force between the drop and the inter- 
face. That force acts mainly on the "barrier ring" of the 
0/W/O film due to heterogeneity of the film thickness 
and the outward flow generated by the collision of the 
drop with the interface. The analytic treatment of the 
problem leads to a probability distribution of rest times 
which can be nicely fitted to the experimental data when- 
ever surfactant is present in the system. Systems without 
added surfactant show a considerable degree of irrepro- 
ducibility. 

The results presented in the previous sections reinforce 
the belief that it is the hydrodynamic resistance and not 
the interaction forces that determine the characteristic 
variation of the coalescence time of the drops as a func- 
tion of their radii. It was shown that the deformation of 
the drops between 100 /xm and 1000 /im, promotes the de- 
velopment of a significant repiilsive barrier which grows 
as a function of the particle radius. Despite this fact, the 
buoyancy force is so large that it suppresses the effect 
of the potentials on the lifetime of the drops. When this 
happens, the behavior found is necessarily determined by 
the mathematical form of the hydrodynamic tensor. 

Up to this point the results correspond to the case 
in which a steady thinning of the intervening liquid be- 
tween a drop and the interface occurs. Therefore they 
cannot explain the large scattering of lifetimes shown by 
the large drops. We explored the possibility that the re- 
ferred scattering could be caused by the contact of the 
capillary waves formed at each oil/ water interface of the 
film. Thus, every time that a drop enter the regions of 
deformation II and III, the program starts to measure 
the lifetime of the film, Tij . Coalescence occurs when 



'^TOTAL 



exp 



(30) 



analytical equation of Vrij [4] for the fastest increase of 
surface oscillations: 



Tvrij = QGn'^j'nho^H^ 



(32) 



Eq. (32) was deduced assuming that only van der 

Waals forces occur between the drops. As a consequence, 
it is generally small and the methodology docs not have 
an appreciable influence on the coalescence time of drops 
of micron size. However, the value of Tvrij depends sen- 
sibly on the thickness of the films, which is large for inter- 
mediate and large drops. It is important to remark that 
in the simulations Tvrij is recalculated at each iteration. 
Therefore, it changes continuously with h. 

Unfortunately, preliminary calculations indicated that 
the methodology described was unsuccessful in reproduc- 
ing the experimental scattering of the lifetimes of the 
drops when the potentials shown in Figures 2 and 3 are 
used in the simulations. It was already shown that the 
buoyancy force is greater than the extensional (and the 
bending) force which occur within region II. Hence, there 
is no repulsive barrier to prevent the rapid thinning of the 
film. When the mechanism of capillary waves is imple- 
mented, coalescence occurs at an even faster rate, but 
this rate is reproducible and does not give rise to a con- 
siderable scattering of lifetimes. 

It is clear that the scattering of lifetimes can only oc- 
cur if a substantial repulsive force overcomes the inertia 
of the particles caused by the buoyancy force. This is 
necessary independently of the mechanism of coalescence 
proposed. Out of all the potential contributions studied, 
it is the steric contribution the one that can be less ac- 
curately determined. The analytic form of the potential 
used in the previous (;alculations is the one formerly em- 
ployed in the estimation of the lifetimes of micron-size 
drops [10]. Thus, it might not be suitable for the simula- 
tion of large deformable drops. Hence, a different expres- 
sion of a steric potential with a harder repulsive barrier 
was tested. This equation results from the adjustment of 
the potential proposed by Alexander-De Gennes [43, 44] 
to the case of truncated spheroids [14]: 

Vst = nr1f{h)+A-KrikBTT''/^Llx (33) 
1.37/ig - Q.2lhf/'^ + 3.20/1"^/'' - 4.36] , 



\hf'^ + -,hr - 1.37 



where hg = h/{2Lg), Lg = Ns, 



ri2 



,1/3 



A^seg is the 



is larger than the width of the film. In the equation above number of segments of the protein (580 aminoacids [31 



Xi{t) = Ran(t) /icrit- 



(31) 



Here, Ran(t) stands for a random number between 
— 1.0 and 1.0. The actual lifetime of a film (r^j ) between 
drops i and j is compared -at each time step- with the 



iseg is the length of each segment (~ 3.0 x 10~ "m) and 
r is the number of molecules per unit area (9.08 x 10^ 
molecules/m^). The value of l^^g was approximated by 

Va^^, where Va is the typical volume of an amionoacid 
residue (K = 57 - 186 A^) [45]. 

It was observed that the lifetime predicted by the cal- 
culations depends considerably on the average length of 
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a protein segment. For Igcg = 5.8 x 10""'^'' m the total 
potential becomes repulsive at short distances and the 
drops do not coalesce even if the mechanism of capillary 
waves is activated. In order to study the influence of the 
new potential on the mechanism of capillary waves, an 
approximate value of Igeg was used {Iseg = 3.0 x 10"^*^ m). 
The value of Iseg was not adjusted to reproduce the order 
of magnitude of the lifetime of the drops. We just substi- 
tuted the steric potential used in the previous simulations 
(section IV A) for the one of Alexander-de Gennes (Eq. 
(33)) to appraise the resulting effect. The steric potential 
used is illustrated in Fig. 13. 



5e+07 




100 



FIG. 13: Steric potential. Solid line: steric potential used of 
the simulations of the Figs. 4-12; dashed line: steric potential 

of Eq. (33) with hag, = 1.85A; dot-dashed line: steric poten- 
tial of Eq. (33) with Zaeg = 3.00^4; dotted line: steric potential 
of Eq. (33) with k^g = 5.80l. 

As expected, it was observed that the lifetime of the 
drops increased with the strength of the steric poten- 
tial. For drops of Vi = 200 /i^m, 400 /xm, and 600 /xm, 
the value of r increased in 15.8, 17.2, and 13.3 seconds, 
respectively. However, the augment of the standard de- 
viation was not sufficient to justify the large dispersion 
of the experimental data. The deviation increased from 
0.0017 to 0.0094 for r, = 200 ^m, from 0.0019 to 0.020 for 
n = 400 /xm, and from 0.0039 to 0.019 for n = 600 /xm. 

In Emulsion Stability Simulations, the interaction po- 
tentials depend on the degree of surfactant adsorption 
to the 0/W interface. Previous calculations [8, 46] con- 
firmed that the coalescence time between two drops is 
inversely proportional to the surfactant surface excess at 
the interface of the drops. Due to the characteristics of 
the experimental set up under consideration, it is very 
likely that drops with different surface excesses might be 
produced. It was already shown [10] that the dispersion 
of the lifetime of small drops is caused by their Brownian 
movement. However, the lifetime of small drops is more 
reproducible than the one of large drops. This coincides 
with the fact that small drops arc produced by means 
of emulsification, while large drops are generated in situ. 



Thus, the protein concentration of the large drops is ex- 
pected to grow as a function of time, starting from the 
moment the drop is formed within the aqueous phase. 

In Ref. [46] we studied the problem of time-dependent 
adsorption and its influence on the coalescence time. Lig- 
gieri et al. [47] demonstrated that several cases of mixed 
adsorption kinetics can be reformulated into a diffusion- 
controlled formalism [48] if the diffusion constant of the 
molecules (Dm) is substituted by an "apparent" diffusion 
constant Dapp- 

Table VII shows the results of additional simulations in 
which the effect of a time dependent adsorption was con- 
sidered. In these calculations, the Alexander-De Gennes 
potential [43, 44] for truncated spheroids [14] was used 
(Eq. (33) with l^^g = 3.00l). The radius of the smah 
drop was kept fixed at r, = 200 ^m, and the value of 
-Dapp was changed in several orders of magnitude. The 
mechanism of coalescence that includes capillary waves 
was activated. Using these conditions, the surface excess 
of proteins at the interface of the drops was evaluated 
from [48-51]: 



D 



app 



1/2 



(34) 



where fapp is the apparent diffusion constant, Cp is the 
protein concentration, and t is the time elapsed from the 
moment that the small drop is released. The time re- 
quired for complete protein coverage of the drop can be 
calculated from the same equation substituting the sur- 
face excess by the inverse of the minimum interfacial area 
of the protein at the oil/water interface: 



AcCp 



(35) 



9.08 X 10^5 



where tc is the critical time, Fmax 
proteins/m^, and Ac = 2(£>app/7r)^/^. 

Notice that when the surface excess of the protein 
increases as a function of time, the value of the in- 
terfacial tension decreases between 7 = 50mN/m and 
7 = 15mN/m. Hence, hini, rf, and the interaction po- 
tentials also change as a function of time. Moreover, the 
value of rf is not fixed within region HI, but also changes. 

From the results shown in Table VII, the following con- 
clusions can be drawn: 

1. If the apparent diffusion constant of the protein is 

high, -Dapp > 2.95 X lO^'' m^/s, the repulsive poten- 
tial between the drop and the interface builds up 
very quickly. Consequently, the coalescence process 
is substantially slowed down, favoring a lifetime of 
the order of 17.5 seconds. 

2. If £>app < 1-0 X 10~^ m^/s, the coalescence of the 

drop and the interface occurs very fast, r w 1.4 s. 
The protein does not arrives in time to delay the 
thinning of the film. 
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3. Intermediate values of Dapp promote a large disper- 
sion of lifetimes. The drops start to feel the re- 
pulsive potential at different film thickness. Hence, 
their coalescence time changes appreciably. More- 
over, the standard deviation of 300 simulations be- 
comes very large, approaching — and surpassing in 
some cases — the order of magnitude of the average 
lifetime. The variation of the deviations indicates 
that for 2.86 x 10"^ mVs < L>app < 2.89 x IQ-^ 
m^/s, slightly different paths might lead to consid- 
erable differences in the lifetime of the drops. These 
are caused by different values of r{t). 



TABLE VII: Variation of the r whit the apparent diffussion 
constant, a represents the standard deviation. 





(mV 


3) 


T (S) 


a(s) 


tc (S) 


1.000 X 10" 


12 


1.45 


7 


50 X 


10" 


-4 


4.90 X 10* 


1.000 X 10" 


10 


1.44 


7 


13 X 


10" 


-4 


4.90 X 10^ 


1.000 


X 


10" 


-9 


1.43 


9 


19 X 


10" 


-4 


49.02 


2.000 


X 


10" 


-9 


1.56 


6 


88 X 


10" 


3 


24.51 


2.250 


X 


10" 


-9 


1.65 


9 


60 X 


10" 


-3 


21.79 


2.500 


X 


10" 


-9 


1.81 


1 


48 X 


10" 


-2 


19.61 


2.750 


X 


10" 


-9 


2.15 


3 


32 X 


10" 


-2 


17.83 


2.800 


X 


10" 


-9 


2.30 


5 


03 X 


10" 


-2 


17.51 


2.850 


X 


10" 


-9 


2.58 


9 


95 X 


10" 


-2 


17.20 


2.860 


X 


10" 


-9 


2.94 




1.91 




17.14 


2.865 


X 


10" 


-9 


4.18 




4.38 




17.11 


2.870 


X 


10" 


-9 


7.02 




6.66 




17.08 


2.875 


X 


10" 


-9 


10.17 




7.34 




17.05 


2.880 


X 


10" 


-9 


13.29 




6.64 




17.02 


2.885 


X 


10" 


-9 


14.95 




5.55 




16.99 


2.890 


X 


10" 


-9 


15.98 




4.46 




16.96 


2.900 


X 


10" 


-9 


17.49 


9 


98 X 


10" 


-2 


16.90 


2.950 


X 


10" 


-9 


17.49 


9 


98 X 


10" 


-2 


16.62 


3.000 


X 


10" 


-9 


17.49 


9 


96 X 


10" 


-2 


16.34 


3.250 


X 


10" 


"9 


17.49 


9 


93 X 


10" 


-2 


15.08 


3.500 


X 


10" 


'9 


17.49 


9 


92 X 


10" 


-2 


14.01 


3.750 


X 


10" 


-9 


17.48 


9 


99 X 


10" 


-2 


13.07 


4.000 


X 


10" 


-9 


17.48 


1 


01 X 


10" 


-1 


12.26 


6.000 


X 


10" 


-9 


17.46 


9 


94 X 


10" 


'2 


8.17 


8.000 


X 


10" 


-9 


17.45 


9 


69 X 


10" 


-2 


6.13 


1.000 


X 


10" 


-8 


17.44 


9 


32 X 


10" 


-2 


4.90 


1.000 


X 


10" 


-7 


17.28 


9 


21 X 


10" 


-2 


4.90 X 10"^ 


1.000 


X 


10" 


-6 


17.28 


9 


21 X 


10" 


-2 


4.90 X 10"^ 


1.000 


X 


10" 


-5 


17.28 


9 


21 X 


10" 


-2 


4.90 X 10"^ 



It appears then that the spreading of the lifetime of 
the large drops could probably be justified in terms of a 
time-dependent protein adsorption. It must be kept in 
mind however, that in our simulations the same value of 
the surface excess is used for the drops and the interface, 
while in the experiment, it is only the approaching drops 
which are likely to be partially covered by the proteins. 



Ghosh and Juvekar [42] published photographs which 
demonstrate that the collision of a millimetric-size (ji = 
1.6 mm) drop with a planar interface resembles the be- 
havior of an under-damped system: "the drop and the 
interface undergo an oscillatory (up-and-down) motion 
before attaining the rest position" . As discussed by these 
authors, elevations of the drop during this motion with 
respect to its final equilibrium position at rest, indicate 
that the movement is similar to that of a rubber ball 
dropped on a stretched membrane. Such behavior is very 
different from the over-damped motion predicted by film- 
thinning models. 

In regard to the above observations, we studied the mo- 
tion of the drops assuming a time-dependent surfactant 
adsorption but deactivating the mechanism of capillary 
waves. It was observed that: 

a. If a very small diffusion constant is used, -Dapp = 

2.9 X 10~^^ m^/s, the film thins smoothly until co- 
alescence occurs. 

b. If a very large diffusion constant is used, Dapp = 2.9 x 

lO"*" m^/s, the film thins monotonically until it 
reaches h = 29.6 nm. After this time {t = 1.19s), 
the drop maintains its average distance from the 
interface fluctuating within 0.30 nm. 

c. For intermediate values of -Dapp like 2.9 x 10"'' m^/s, 

a drop of = 200 fiia reaches much smaller separa- 
tions h = 11.8 nm {t = 1.6 s) before the repulsive 
barrier builds up. Next, the drop is pushed out- 
ward by the potential until it reaches h ~ 29.7 nm. 
Following it oscillates around this distance. 

d. When the mechanism of capillary waves is activated 

for the case of intermediate apparent diffusion con- 
stants (i.e. -Dapp = 2.9 x 10~^ m^/s), the coa- 
lescence occurs during the outward motion of the 
drop. It does not happen at its closest separation 
from the surface. 

From the above observations it appears likely that the 
bouncing of a drop at the oil/ water interface is caused 
by the sudden appearance of a repulsive force at a short 
distance of approach. That repulsive force might occur 
as a result of a time-dependent adsorption, but it might 
also happen through the accumulation of surfactant at 
the so-called "barrier ring" as proposed by Ghosh and 
Juvekar [42]. 

V. CONCLUSION 

The lifetime of drops pressed by buoyancy against a 
planar interface can be reproduced using Emulsion Sta- 
bility Simulations. Our calculations support previous 
theoretical evidence that suggests that the particular 
shape of the curve of r vs. is basically determined by 
the drainage of the intervening film between the drops 
and the interface. However, depending on the particle 
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size, different expressions of the initial distance of de- 
formation and the maximum film radius of the drops 
are necessary. In the case of large drops the appropri- 
ate equations results from the consideration of the buoy- 
ancy force experienced by the drops. In the case of small 
non-deformable drops the (attractive) potential of inter- 
action plays a significant role. The intermediate range 
of particle sizes between these two extremes results from 
a balance between interaction forces and the buoyancy 
force. Still, the calculations support the occurrence of 
deformable drops between 10 /xm and 100 /xm of particle 
radii. However, further study of this size range is neces- 
sary. 

While the average value of the coalescence time can be 
justified by the mechanism of film thinning, the disper- 



sion of the lifetime of large drops cannot be. As shown 
by the simulations, a possible explanation of this phe- 
nomenon comprises a combination of a substantial re- 
pulsive barrier for coalescence, the occurrence of capillary 
waves and a time-dependent surfactant adsorption. 
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